Timing Characteristics of the Hard X-ray Emission from Bright 
X-ray Pulsars Based on INTEGRAL Data 



© 2009 A.A.Lutovinov 1 *, S.S.Tsygankov 21 

1 Space Research Institute, Russian Academy of Sciences, Profsoyuznaya ul. 84/32, 

Moscow 117997, Russia 
2 Max- Planck- Institut fur Astrophysik, Karl- S chwarzs child- Str. 1, Postfach 1317, D-85741 

G arching, Germany 

Received June 17, 2008 

We review the results of a timing analysis of the observations for ten bright X-ray pulsars 
(with fluxes > 100 mCrab in the 20-100 keV energy band) that fell within the INTEGRAL 
field of view from 2003 to 2007. The dependence of the pulse profile on the energy and intrin- 
sic source luminosity has been investigated; particular attention has been paid to searching 
for changes in the pulse profile near the cyclotron frequency. The dependence of the pulsed 
fraction for X-ray pulsars on their luminosity and energy band has been studied in detail for 
the first time. 
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INTRODUCTION 



According to the universally accepted theory of accretion onto rapidly rotating neutron 
stars with a strong magnetic field, the infalling matter from the companion star is braked on 
the Alfven surface by the magnetic pressure, becomes frozen into the magnetic field, and falls 
along the field lines into narrow ringlike regions at the magnetic poles of the neutron star, 
producing accretion columns and releasing its gravitational energy in the X-ray wavelength 
range. In that the pulsar rotates about its axis, an observer will register pulses on the 

light curve differing in shape and features, depending on the specific physical and geometrical 
conditions both near the neutron star surface and in the path of signal propagation. The 
pulsed fraction will depend on the configuration of the emitting regions, the dipole position 
relative to the observer, the energy, etc. 

As was pointed out in early papers on this subject (see, e.g., Wang and Welter 1981; 
White et al. 1983), the observed pulse profiles differ significantly for different sources, have 
a wide variety of shapes, depending on the energy and source luminosity, and, in some cases, 
can shift by as much as 180° in phase as the energy band changes. In addition, they are 
also variable on the scale of a single pulse (see, e.g., Frontera et al. 1985; Tsygankov et al. 
2007). Nevertheless, Bulik et al. (2003) made an attempt at a "broad" classification of pulse 
profiles, according to which each pulsar can be assigned to the subclass of either single- or 
double-peaked profiles. This classification is based on the fact that an observer sees either 
one or both magnetic poles, respectively, as the neutron star rotates. This difference may not 
be so distinct at low energies, but single- and double-peaked profiles can be clearly separated 
above ~ 10 keV, where the effect of absorption may be neglected, or, in the case of pulsars 
with a cyclotron feature, above the cyclotron energy. 

A strong magnetic field (10 11 — 10 13 G) near the emitting regions on the neutron star 
surface is the source of peculiar features in the observed properties of X-ray pulsars. In 
particular, the position of the cyclotron resonance scattering feature in the pulsar spectrum 
is a direct source of information about the magnetic field strength. In this way, the magnetic 
field on the X-ray pulsar Her X-l was measured for the first time (Truemper et al. 1978). 
In addition, the properties of the accreting plasma can change sharply at the cyclotron 
frequency, which can be manifested as a change in the beam pattern of emission (Meszaros 
and Nage 1985). The corresponding changes in the pulse profile near the cyclotron energy 
were observed in some of the sources (see, e.g., Tsygankov et al. (2006) and references 
therein) . 

Basko and Sunyaev (1975, 1976) showed that the beam pattern of emission from an 
accretion column depends on the presence of a shock in it. Thus, for example, a shock 
in which the infalling matter is braked emerges above the neutron star surface at a high 
luminosity (> 10 37 erg/s). The emitting plasma is accumulated in the zone under the shock 
and the emission emerges mainly through the side surfaces of a accretion column with a 
fan beam pattern. At lower luminosities, the matter can be braked near the neutron star 
surface and a pencil beam pattern of emission will be more likely due to the effects of a 
strong magnetic field. At intermediate luminosities, the beam pattern will be a combination 
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of fan and pencil beams. The beam pattern and, accordingly, the pulse profile can change 
with energy band, since the hotter layers of matter are closer to the neutron star surface, 
where the emission formation conditions can differ from those in higher layers (Basko and 
Sunyaev 1976). 

Obviously, changes in the beam pattern and in the emission formation conditions and 
geometry will lead to a dependence of the pulsed fraction on both source luminosity and en- 
ergy. This was recognized and measured by several authors back in the 1980- 1990s. However, 
only with the advent of the modern RXTE and INTEGRAL observatories, with high time 
and energy resolutions (especially in the hard energy band, where the observed emission 
is unaffected by photoabsorption and depends only on the system geometry and physical 
conditions in its formation region), did systematic studies become possible. In particular, 
using RXTE data, Tsygankov et al. (2007) showed that the pulsed fraction for 4U01 15+63 
decreases with increasing source luminosity and increases with energy. The increase is not 
monotonic but has features near the harmonics of the cyclotron absorption line: its excess 
above the general trend is clearly observed near the fundamental harmonic; similar features 
are also observed near the higher harmonics, but they can be slightly shifted relative to 
the center of the cyclotron line harmonic. An increase in the pulsed fraction with energy 
was also reported for several other pulsars: GX 1+4 (Ferrigno et al. 2007), OAO 1657-415 
(Barnstedt et al. 2008), and EXO 2030+375 (Klochkov et al. 2008a). 

This is the next paper in our series of papers (Tsygankov et al. 2007; Lutovinov and 
Tsygankov 2008) devoted to the investigation of the pulse profiles and pulsed fractions for 
X-ray pulsars. Using the methods proposed and described by Tsygankov et al. (2007), we 
investigated the ten (4U 0115+63, V 0332+53, A 0535+262, Vela X-l, Cen X-3, GX 301-2, 
OAO 1657-415, Her X-l, GX 1+4, EXO 2030+375) brightest (with fluxes > 100 mCrab) 
X-ray pulsars based on INTEGRAL data in the hard (> 20 keV) energy band. Here, we 
analyze in detail the dependence of the pulsed fraction for X-ray pulsars on their luminosity 
and energy band and the effect of resonance absorption near the cyclotron line on it for the 
first time. In addition, we obtained pulse profiles for the above pulsars in several energy 
bands, relative intensity maps of pulse profiles, and dependences of the pulsed fraction 
(PF) on energy for all of the possible observed states and luminosities and compiled the 
corresponding catalog. Since this is a fairly large data set and in order not to overload 
the paper, we provide here only the averaged (if the profile did not change significantly with 
luminosity) or most typical pulse profiles, intensity maps, and PF dependences on the energy. 
The complete catalog is freely accessible at http://hea.iki.rssi.ru/integral/pulsars and can 
be used to construct and test radiation models for X-ray pulsars. 

OBSERVATIONS AND DATA ANALYSIS 

Data from the IBIS telescope (the ISGRI detector; Lebrun et al. 2003) of the INTEGRAL 
observatory (Winkler et al. 2003) were used to investigate the pulse profiles in the hard X-ray 
energy band. 
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Table 1. Bright X-ray pulsars from INTEGRAL data 



Pulsar 


Observation, 


Middle 


Exposure, 


Flux a , 


Luminosity^* , 


name 


orbits 


of observations, MJD 


ks 


xlO — 9 erg s — 1 cm -2 


10 37 erg s — 1 


l 


2 


3 


4 


5 


6 


4U 0115+63 


238 


53273.6 


97.0 


12.9+2^ 

—0.2 


7.6+°-? 

1 "—0.1 




272 


53376 5 


23.7 


78 4+2i 


46 1+2-1 




273 


53379.3 


42.5 


74.0+2? 

— U. / 


43.5+^ 




274 


53380.3 


15.9 


70.2+°^ 

—0.9 


41.2+21 

—0.5 


V 0332+53 


278 


53394.2 


41.9 


45.3+°' 

— O.o 


26.6+21 

U. 2 




284 


53411.0 


107.7 


15.0+°? 

—0.1 


8.8+°+ 

" — 0.1 




285 


53413.1 


11.6 


14.8 + °^ 

— 0.3 


8-7+21 

— U.z 




286 


53416.3 


16.2 


9.9l|? 


5.81^1 




287-288 


53420.6 


27.6 


K q+4-6 

—0.8 


3.4+^1 

" — 0.4 




352 flow) 

'.J " . / \ Ivy vv / 


53613 7 


34.1 


15 6+2! 
±u - u — 0.1 


1 3+2'} 

0.1 


A 0535+262 


352 (medium) 


53614.5 


51.0 


18.8+°! 

—0.1 


1.5+21 

—0.1 




352 (high) 


53615.4 


42.0 


20.9+2-1 

— O.O 


1.71°-} 




58 


52734.3 


7.2 


10.7+^ 

— U.O 


o.5i+°-°J 

— O.Uo 




81 


52804.1 


60.0 


5 63+2^ 


o 27+2-2! 

w 1 —0.01 




82 


52807.5 


106.5 


q qq+o.16 
— u.zy 


0.45+2-2} 

— U.U1 




83 


52809.5 


150.8 


9 23+ - 12 

— U. lo 


0.44+°-°! 

— U.U1 




84 


52812.7 


61.5 


4.53+2-Ji 

—0.18 


0.22+2-2! 

—0.01 




85 


52816 7 


79.9 


8 07+S-15 

u,w 1 —0.24 


o 39+2-2! 




86 


52819 1 


148.1 


8 93+2 of 


o 43+2-2! 

W - 4: ' J -0.01 




87 


52822.3 


57.8 


4 82+2^ 

^•"^-0.23 


o 23+2-2! 

"•■ 4 " J — 0.01 




88 


52825.0 


88.6 


12 60+2'os 
- L ^ i,uw — 0.28 


o 61+2-2! 




137 


52971.4 


111.1 


11 68+2^ 


o 56+2-2! 

0.02 




138-140 


52978 5 


248.1 


8 64+2™ 

o.<j"±_ q 28 


o 42+2-2! 

kJ-^^ — 0.01 




146 


52997.5 


5.2 


18.6+9-I 

u 1.6 


0.9+2-! 

" — 0.1 




149 


53006.4 


5.5 


20.0+J-jj 


i.o+g-1 

— U.o 




154 


53021 3 


5.9 


14 6+2-f 


o 7+2-! 

u - ' -0.1 




157 


53030 3 


8.7 


10 i+?i 

±U -- L — 1.8 


o 5+2-! 


Vela X-l 


161 


53042 2 


6.8 


7 64+!-" 

' • u ^ t — 1.78 


o 37+2-2^ 

U -' J 1 -0.09 




171 


53071 6 


5.0 


11 6+2« 

— 0.6 


55+2 ™ 

v.u<j_q Q 3 




186 


53119.2 


5.5 


— U.o 


0.58+°-°4 

— U.U4 




203 


53168.5 


9.8 




55+2^ 

" • KJKJ — 0.13 




217 


53209.7 


7.1 


12.5+° : ? 
"—l.i 


0.60+2^ 

—0.05 




250 


53308 5 


8.0 


8 98 + no? 


43+2-2^ 




256 


53326 6 


8.4 


12 2+?^ 

— 1 9 


58+2no 




262 


53344.5 


7.6 


19.4+° : ? 


0.93±g-Sg 




267 


53361.4 


28.4 


10.718J 


o-5i±8:g 




275 


53383.6 


10.4 


7.561^ 


0.361HI 




301 


53460.9 


4.9 


i o q+0.5 
io.o_ g 


0.64±g;S 




373-383 


53693.1 


783.0 


9.661°;!° 


0.46+°;°! 
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1 


2 


3 


4 


5 


6 




398 


53751.1 


4.7 


23 7 +0 n 3 A 

Z " J - 1 -0.4 


1 1 0+0.01 

1 - ld -0.02 




433-440 


53865.9 


749.5 


8.58 + °-|° 

— U.lo 


0.41+H1 

— U.01 




445 


53891.6 


24.7 


6 

" ,UJ -0.34 




"•^-0.02 




446 


53894.8 


7.0 


6.85^ 


0.33+g-gi 

— U.Uo 




192 


53136 1 


94.5 


10 6+^ 

lu,u -0 .3 


°-- L — 0.2 




194-197 


53146 8 

.J '.J _L -L V7 • V7 


349.3 


5 05+2-^ 

u,uu — 0.13 


3 88 +0 /! n 
d,uu — 0.10 




198-200 


53157.0 


264.0 


5.43 + °-°I 

—0.04 


—0.03 


Cen X-3 


201-205 


53165.5 


175.8 


4 16 +0 / n7 

^" XVJ — 0.07 


3 20+S-Sf 

<j.^«j_ 0.05 




325 


53533 9 


44.7 


<j.u-*_ o.o4 


9 7 9+0.09 
z ' ' -0.03 




327 


53539 9 


80.1 


2 04 +0 ,- nQ 


x,d 1 —0.07 




330 


53549.0 


62.6 


5.37+gi? 

— U. 1 1 


4 19+0.09 

— U.Uo 




46 


52700 5 


4.0 


11 7 + °^ 

±± - ' —0.2 


3 

d - ^-0.06 




76 


52789 5 


44.8 


3 77 + n ?I 
°- ' ' -0.11 


-"-•^ 1 -0.04 




77 


52791.1 


29.5 


i,uu -0 .05 


68 +0 ,-^ 

u.uo_o.02 




78 


52796 

.7 1 '7 V7 • V 7 


24.6 


4 19 +0 ,no 

^ - x J — 0.09 


- 1 - •^ t - L — 0.03 




91 


528343.0 


34.2 


2.65 + °-li 

— U. 14 


0.89+2$ 

—U.Uo 




150 


53009 8 


33.4 


^- 1 -0.07 


1 59+n no 
±.o^_o 02 




176 


53087 4 


103.9 


2 oo+S-Sl 

z,.uu_ 04 


67+°,^ 

u - u ' -0.01 




213 


53198 


3.3 


2 Q1+ - 02 
^ - ^ — 0.03 


98 + S-Si 

u.^o_q 01 




267 


53361.1 


3.6 


4 7O+0.08 
—0.09 


1.62 + °$ 

—0.03 


GX 301-2 


268 


53362 5 


4.9 


8 09+°-?? 


9 70+O.O3 
z - ' d -0.05 




283 


53407.4 


3.2 


12. 5_ 02 


4 22+°$ 

^•^^—0.06 




322 


53524 9 

.J '.J .J J- • '.7 


112.4 


oo+O.Ol 
^.<j>u_ q 03 


—0.01 




323 


53528 2 


104.4 


7 07+°-°i 

' - u ' -0.05 


2 38+°,-°,i 

z,.oo_ 02 




325 


53533.4 


31.4 


0.00 n ri7 
— U.U < 


l-91 + o'8l 

— u.uz 




326 


53536 8 


73.3 


2 70 +0 ,ni 

1 u -0.04 


91+°,m 
* — 0.01 




330 


53549.4 


63.6 


3.71+2-Si 

1 -"-—0.05 


1 9 f-+0.01 
"—0.02 




407 


53778.4 


3.1 


11 9 + °J 

± - L - E '— 0.3 


4 02 +0 r? 5 
^• Ui -o .11 




444 


53888.8 


3.1 


3.60+°-?^ 

" uu -0.13 


1 91+0.03 

—0.05 




463 


53946.1 


3.1 


10 6 +0 / 9 


3 59 + n -o7 




517 


54108.5 


8.6 


6.01+0 1 - ? 

— U.ll 


2.02+^4 

— U.U4 




36-38 


52671 8 

.J V7 1 _J_ - v 7 


142.1 


3 26 + ^-? 7 

<j.^vj_0.17 


1 60+°,°,« 

x ,VJW — 0.08 




46-47 


52701.1 


145.4 


O 7 q+0.07 

°- 1 d -0.13 


1 OO+0.03 

i.Od_o 06 




50-52 


52714.6 


169.4 


3 69 +0 ,-?l* 


1 01 +0.03 

y-, 01-. 

x - ox -0.08 




60-61 


52742.8 


66.3 


4 13+°-?^ 

±,J -0.17 


2 03 + °-°f 

z,.w<j_Q 08 




103 


52870.0 


110.6 


5.03^;^ 


2.471^2 




105-110 


52883.3 


151.9 


pro+0.09 
O.OO_ .20 


1 7fi +0.04 
i -'°-0.10 




116 


52908.4 


48.6 


4 qo+o- 08 

^• yz -0.24 


9 49+0.04 
z -^ z -0.12 




119 


52916.7 


53.1 


oq+0.08 
O.OO_ .27 


1 oo+0. 04 

i-oo-o.^ 




175 


53085.2 


36.4 


1 qpr+0.06 


66+ 03 



1 


2 


3 


4 


5 


6 




224 


53232.5 


75.6 


4 QO+0.07 


9 1 9+0.03 
z - ±z -0.06 




233 


53256.9 


17.9 


1 n 1+0.06 
i - ui -0.12 


o.5o±8:8| 




283-284 


53409.9 


72.2 


9 99+O.O6 
z - zz -0.15 


i.09i8:8 3 7 




285-289 


53420.4 


231.7 


q ci +0.08 
°- oi -0.16 


1 7Q+0.04 
L - ' °-0.08 


OAO 1657-415 


290-294 


53434.1 


230.0 


9 £9+0.05 
z - uz -0.11 


1 oq+0.03 




295-299 


53450.5 


182.9 


q oq+0.08 
O.OO_ 21 


1 OO+0.04 
-L-Sci-O.lO 




303-307 


53474.3 


168.3 


o oq+0.06 
°- oy -0.20 


1 67+ - 03 

L - v ' -0.10 




345 


53594.9 


24.5 


i i 7+0.03 
i - i '-0.65 


U -°°-0.32 




347 


53599.7 


81.2 


9 a q+0.06 
z -^°-0.09 


1 1 q+o.03 

i - ±y -0.04 




350 


53608.7 


91.4 


i 91+O.O3 
1 - z,1 -0.06 


o.6o+°:°I 




364 


53650.1 


83.6 


6.1918:3 


3-04±8:?l 




409 


53785.1 


69.1 


A 97+O.O8 
' -0.26 


2 01+ - 04 

z - u± -0.13 




412 


53794.8 


44.7 


9 70+0. 08 
z - ' °-0.14 


1 04+O.O4 
i -°^-0.07 




472-473 


53974.8 


56.4 


°- uo -0.30 


9 q7+ - 07 

z - y ' -0.15 




525 


54131.4 


91.6 


5 96 +0 - 12 
o.yu_ 24 


9 Qq+0.06 
z,.yo_ 12 




531 


54149.5 


79.7 


oo+0. 08 
°- oo -0.13 


1 66+ - 04 
i - uu -0.06 




536 


54164.6 


81.7 


2 Q1+ - 07 

z - yl -0.12 


1 43+0.03 

i -^°-0.06 




540 


54176.1 


49.4 


a 7+O.O6 
-0.17 


1 71+O.O3 
L - ' i -0.08 




339 


53575.8 


132.3 


«c;+0-02 

°- ud -0.09 


1 co+0.01 
l.UO_ 4 


Her X-l 


340 


53577.8 


19.5 


5.0918:" 


9 on+0.01 
z " z ' u -0.10 




341 


53582.4 


40.4 


4.iol8:S 


1 77+O.07 
L ' ' -0.09 




342 


53584.2 


67.9 


1 7c;+o.ud 
L - ' °-0.10 


u - ' °-0.04 




46 


52699.4 


100.0 


q crcr+0.01 


1 cq+0.01 
± - JO -0.05 




60-64 


52746.9 


149.3 


fiS+ 01 

u -°°-0.03 


o.27i8:8i 




119-122 


52923.1 


283.3 


1 qo+0.01 
± - oo -0.04 


5Q+ - 01 
u - ja -0.02 




164-165 


53053.8 


207.2 


1 oq+0.02 
± - oo -0.06 


7Q+ - 01 
u - ' y -0.03 




168 


53064.1 


74.0 


O 70+O.OI 


1 1 9+O.OI 
± - ±z -0.05 




172-173 


53078.5 


152.4 


1 90+O.O2 


u -°°-0.02 




175 


53084.9 


38.9 


1 91+O.OI 
i - zi -0.07 


u - oz -0.03 




178-179 


53095.3 


231.4 


1 04+O.OI 


o.8oi8:8 2 


GX 1+4 


181 


53103.7 


95.4 


9 90+O.O2 
z - zo -0.06 


0.9918:83 




183 


53108.5 


48.3 


1 01 +0.01 
± - oi -0.07 


p;7+ - 01 

v -° ' -0.03 




185 


53114.8 


29.9 


4 02 +4 - 02 

^• uz -0.12 


1 74+ 1 - 74 

l - '^-0.05 




225-234 


53248.5 


691.2 


2 84 +0 - 01 

z -°^-0.05 


1 9 q+0.01 

1 - zo -0.02 




235-237 


53267.3 


291.7 


4^+0.01 

°-^°-0.05 


1 4 q+o.oi 

i -^ y -0.02 




240-241 


53281.1 


250.8 


O /IQ+0.01 

<J- 4y -0.06 


1 r. 1+0.01 
L - ol -0.02 




242-246 


53291.2 


238.6 


1 fi o+0.01 


7n+ - 01 

u - ' u -0.03 




297-299 


53453.7 


158.1 


n fji +0.01 
u - ui -0.03 


u - zu -0.01 




303-305 


53471.3 


181.4 


i -°^-0.04 


0.7118:82 




306-308 


53479.2 


277.9 


i - yi -0.06 


0.8318:8s 




475-477 


53985.2 


211.5 


2 fiS +0 - 01 

z - uo -0.03 


i-i4l8:81 
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1 


2 


3 


4 


5 


6 




67 


52761.8 


64.4 


o p-9+0.03 
z - uz -0.04 


i p-9+0.02 
1 -°^-0.03 




80 


52801.2 


129.2 


1 yo+0.03 
L - ' °-0.04 


1 ns+ 002 

i - uo -0.02 




159-160 


53040.9 


249.4 


2 40 +0 - 01 

z -^ u -0.04 






190 


53130.2 


122.8 


o on+0-02 
^•° U -0.05 


9 90+O.OI 
z - zo -0.03 


EXO 2030+375 


218 


53214.2 


99.0 


i -i O+0.03 
i - io -0.04 


o-68±8:S 




251-255 


53319.2 


339.4 


o nO+0-02 
°- uz -0.04 


1 01+0.01 
J-^ 1 - 0.02 




437-438 


53870.3 


126.4 


i-97±8:S 


1 1 0+0.01 

i - i8 -0.03 




452-453 


53914.9 


16.5 


2i-7±8:i 


i3.ol8:I 




462 


53942.9 


70.1 


29.31™ 


i7.eir 7 




470 


53967.6 


140.1 


28.3+°i 


17.0183 




486 


54015.2 


133.6 


4.29^ 


9 fil+0.18 
Z - O1 -0.24 



In the 3-100 keV energy band. 

In the 3-100 keV energy band assuming the distance to 
the source to be known (see table 2). 
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The IBIS data were processed for a timing analysis using software developed and main- 
tained in the National Astrophysical Institute in Palermo, Italy (http:/ /www.pa.iasf.cnr.it/~ 
ferrigno/INTEGRALsoftware.html); a description of the data processing technique can be 
found in Segreto and Ferrigno (2007) and Mineo et al. (2006). Since one of our main ob- 
jectives was to search for a correlation between timing and spectral characteristics of the 
emission, we also performed a comprehensive spectral analysis (in particular, we searched 
for and determined the parameters of cyclotron features) for the entire set of sources under 
study using the same software as that applied for our timing analysis. In accordance with the 
user's guide for the INTEGRAL data processing software, the corresponding component at 
a 2% level was added in the XSPEC package to take into account the systematic uncertainty 
in our spectral analysis. 

A log of pulsar observations is presented in Table 1. Its columns give (1) the source 
names, (2) INTEGRAL observatory orbit numbers, (3) dates of observations (in MJD), (4) 
exposure times, (5) source fluxes, and (6) source luminosities (calculated by assuming the 
distance to the source to be known) in the 3 — 100 keV energy band. The published sources 
of distances to the objects under study are presented in Table 2. 



Table 2 



Source name 


Distance, kpc 


Reference 


4U 0115+63 


7 


Negueruela and Okazaki (2001) 


V 0332+53 


7 


Negueruela et al. (1999) 


A0535+262 


2.6 


Janot-Pacheco et al. (1987) 


Vela X-l 


2 


Sadakane et al. (1985) 


Cen X-3 


8 


Krzeminski (1974) 


GX 301-2 


5.3 


Kaper et al. (1995) 


OAO 1657-415 


6.4 


Chakrabarty et al. (2002) 


Her X-l 


6 


Howarth and Wilson (1983) 


GX 1+4 


6 


Chakrabarty and Roche (1997) 


EXO 2030+375 


7.1 


Wilson et al. (2002) 



It should be noted that since different authors used different models (Gaussian or 
Lorentzian profiles) in fitting the cyclotron line, there is a spread in its energies for the 
same sources and their states in the literature. Here, we used a Lorentzian profile in the 
form 



CYCL(E) = exp( ™l^J2* )' (1) 
v {h - h cyd ) + a cyd / 

where E cyc \ is the line center, T cyc \ is the line depth, and o cyc \ is the line width (Mihara et 
al. 1990). In our spectral analysis, we added this component to the "standard" fit for X-ray 
pulsars (a power law with a high-energy cutoff; White et al. 1983). 

The time dependences of the observed source flux in the 20-60 keV energy band and the 
pulse period are shown in Fig. 1. The pulsar 4U 0115+63 was observed by INTEGRAL 
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Time, MJD 



Fig. 1: Light curves for the pulsars under study (in the 20-60 keV energy band; 
crosses) and time dependences of their periods (open circles). The error in the period 
was determined by the method described by Filippova et al. (2004) and corresponds 
to la. 



for a fairly short time interval in September 2004, in which its intensity and pulse period 
underwent no significant changes; for the pulsar Her X-l, the corresponding periods and 
fluxes were obtained by Klochkov et al. (2008b). 
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INTEGRAL/1BIS/ISGRI Crab 




Pulse phase ^ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

30 40 

Energy, keV 



Fig. 2: (a) Pulse profile for the pulsar NP0532 in the Crab Nebula from INTEGRAL 
data in the 20-50 keV energy band, (b) Relative intensity map in "energy — pulse 
phase" coordinates (for more detail, see the text). 

Our main results are the pulse profiles in wide energy channels (20-30, 30-40, 40-50, 50-70, 
70-100 keV), the relative intensity maps in "energy - pulse phase" coordinates from IBIS 
data, and the dependencies of PF (calculated as PF = , m "' , where I max and l m i n are 

J-max\ J-min 

the source intensities at the maximum and minimum of the pulse profile) on energy and 
luminosity. The intensity maps were obtained by convolving the pulsar light curve corrected 
for the background emission in narrow energy channels (the width was chosen to provide an 
optimal significance and was about 5 keV) whose centroid shifted from channel to channel 
by 2-3 keV. Each profile was constructed in units relative to the mean count rate in a given 
channel. We give the resulting map normalized to unity (all intensities were divided by the 
maximum value over the entire map). Such a representation reflects well the evolution of the 
profile structure (the change in the relative contribution of different peaks etc.). For each 
source, the results mentioned above are either averaged if the pulse profile did not change 
significantly with luminosity or most typical for different states (the source state is given 
in the corresponding captions to the figures). The dotted lines in the corresponding figures 
indicate the positions of the harmonics of the cyclotron absorption line (if it is observed in 
the spectrum). Note that the background is difficult to properly subtract at high energies 
(> 50 — 60 keV, depending on the source spectrum and intensity) and, therefore, the formal 
values of PF can exceed 100% in some cases due to insufficient statistics. In current paper 
all relative intensity maps are presented in a black-and-white style; more brighter colors 
correspond to regions of higher intensities. The color maps are available in the electronic 
publication. 

To illustrate the operation of the method used and to estimate the accuracy and validity 
of our timing analysis, we performed the corresponding calculations for the pulsar NP0532 
in the Crab Nebula. This source was observed by INTEGRAL once in several months to 
calibrate the instruments. Figure 2 shows the pulse profile in the wide 20-50 keV energy 
channel and its relative intensity map. The profile has a characteristic double-peaked shape 
with a slight increase in the relative intensity of the second (relative to the zero phase we 
chose) peak with energy, which is clearly seen from the intensity map (Fig. 2b). 
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Fig. 3: (a) Time dependence of the pulse period for NP0532 in the Crab Nebula 
from INTEGRAL data (filled circles); the solid line represents an extrapolation of the 
Jodrell Bank radio measurements, (b) The corresponding deviations of the INTE- 
GRAL measurements from the extrapolation (filled circles) and the radio measure- 
ments (dashed line). In both ranges, a glitch is clearly seen at MJD 53070. 



In Fig. 3a, the filled circles indicate the time dependence of the intrinsic pulsar 
period that was measured and reduced to the Solar system barycenter by the meth- 
ods used here; the solid line represents an extrapolation of the pulsar period measure- 
ments in the radio frequency band by the Jodrell Bank Observatory in the form P = 
0.033522231 + 0.00000003633565 x (T — 52014) where T is the time in MJD (Lyne et al., 
1993; http://www.jb.man.ac.uk/pulsar/crab.html). Figure 3b shows the deviations of the 
periods measured by INTEGRAL from the extrapolation. The dashed line indicates the de- 
viations of the actual period measurements by the Jodrell Bank Observatory from the same 
extrapolation. We see excellent agreement between the radio and gamma-ray measurements. 
Note also the detection of a pulsar glitch by both observatories at MJD ~ 53070. 

Comparison of the data obtained suggests that the pulse period determined for the pulsar 
NP0532 by INTEGRAL coincides with a high accuracy (a typical difference ~ 1 x 10 -9 s) 
with that measured in the radio frequency band. 
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RESULTS 

We studied X-ray pulsars in binary systems of various classes; however, all of these sources 
are united by the fact that a relatively high (> 100 mCrab) flux is observed from them in the 
hard X-ray energy band and one or more harmonics of the cyclotron absorption resonance 
line are present in the spectra of most of them (see, e.g., Coburn et al. 2002; Filippova et al. 
2005). Since our analysis required high statistics, the chosen sources either had long exposure 
times or were investigated during their outburst activity. Below, we briefly present the results 
obtained for each X-ray pulsar; our attention was focused on analyzing the dependences of 
PF on the energy and source luminosity. 

4U 0115+63 

The transient X-ray pulsar 4U 0115+63 was observed by INTEGRAL during an intense 
outburst in September 2004 at the decay phase, when the source luminosity was ~ 7.5 x 10 37 
erg/s. The pulse profile characteristic of such a luminosity level in hard energy channels, 
along with the relative intensity map and the energy dependence of PF, are shown in Fig. 
4. Note that the pulse profile properties and the PF variation over a wide range of pulsar 
luminosities, energies, and time were studied in detail by Tsygankov et al. (2007) based on 
the RXTE data obtained during the 1999 and 2004 outbursts. In particular, these authors 
found that a "wavelike" behavior of the profile with changing energy is observed over the 
entire luminosity range of the pulsar 4U 0115+63. The effect lies in the fact that each line of 
the level of equal intensity of the main peak does not lie at one phase for different energies 
but slightly shifts alternatively in one and the other directions. Interestingly, such a phase 
variability of the main peak is repeatable in energy and its period roughly coincides with 
the energy separation between the harmonics of the cyclotron absorption line in the pulsar 
spectrum (Fig. 4). The detection of a strong dependence of PF on energy and luminosity, 
according to which PF decreases with increasing luminosity and increases with energy, is no 
less important and interesting. In the latter case, this increase is not uniform but has local 
maxima near the cyclotron line harmonics (see, e.g., Figs. 3 and 4 from Tsygankov et al. 
2007). Such a behavior of PF is confirmed by INTEGRAL data (Fig. 4). 

Note that 4U 0115+63 was the first X-ray pulsar for which such peculiarities of PF were 
detected. Subsequently, similar dependences were found for several other X-ray pulsars with 
cyclotron lines in their spectra (see below). 

V 0332+53 

Like 4U 0115+63, the X-ray pulsar V 0332+53 was observed during an intense outburst 
in a wide range of fluxes from ~ 70 to ~ 700 mCrab (20-60 keV energy band), which 
corresponds to a spread in luminosities from ~ 4 x 10 37 to ~ 4 x 10 38 erg/s (assuming 
the distance to the pulsar to be d ~ 7 kpc). Three harmonics of the cyclotron absorption 
line were identified in its spectrum (Coburn et al. 2005); the position of the fundamental 
harmonic is not constant with luminosity (Mihara et al. 1998). Tsygankov et al. (2006) 
were the first to show that the cyclotron energy in the source spectrum decreases linearly 
with increasing luminosity during outbursts. Dramatic changes in the pulse profile with 
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luminosity and energy, especially near the cyclotron feature, were also found. In this study, 
we obtained a set of pulse profiles and relative intensity maps for the entire INTEGRAL data 
set (see http:/ /hea.iki. rssi.ru/integral/pulsars). Here, as an illustration, we present only the 
results obtained for high (orbit 273, Fig. 5) and low (orbit 284, Fig. 6) states, where changes 
in the profile with luminosity are clearly seen. 

Despite the above dramatic changes in the pulse profile, the energy dependences of PF 
are very similar in both general increase with energy and local features (maxima) 

near the cyclotron absorption line (Figs. 5 and 6). Interestingly, their positions are almost 
independent of the source luminosity and are shifted toward the higher energies with the 
center at ~ 32 keV. 

AO 53 5 +262 

The X-ray pulsar A0535+262 was observed by INTEGRAL near the maximum of its 
outburst in August- September 2005. Although the total exposure time was not very long, 
we arbitrarily divided it into three parts, depending on the source luminosity (see Table 1) 
and searched for pulsations and constructed the pulse profiles in each of them. The derived 
pulse profile is a fairly wide (~ 0.8 period) peak with two maxima and an indistinct dip 
in intensity between them (Fig. 7). The width of the dip increases with energy, while the 
amplitude of the second maximum decreases, as can be clearly seen from the relative intensity 
map(Fig. 7). In softer channels, the pulse profile is considerably more complex (Caballero 
et al. 2008). 

The source spectrum exhibits two harmonics of the cyclotron absorption line (~ 45 and 
~ 100 keV); the position of the fundamental harmonic measured using data from the Suzaku 
observatory remains almost constant in a wide range of luminosities (Terada et al. 2006). 
However, Caballero et al. (2008) pointed out a possible change in the position of the cyclotron 
line from ~ 50 keV during the precursor to ~ 46 keV during the main outburst. It is seen 
from the intensity map that the relative contribution of peaks in the profile is changed when 
crossing the cyclotron energy, that is typical for other pulsars from current review. 

The source PF increases with energy; no distinct features are observed near the cyclotron 
line harmonics, which is most likely because the line depth is not so high and the line is 
fairly wide. 

A change in the pulse period of the source with its intensity during the outburst is clearly 
seen from Fig. 1. Such a spin- up of the neutron star can be roughly described in terms of 
the theory for angular momentum transport from the disk to the neutron star (Ghosh and 
Lamb 1979): 

)L 6 / 7 = ^ 7 n(uj s )(4ird 2 Ff 7 , (2) 

where fi is the magnetic moment of the neutron star with magnetic field B and radius 
R, n{oj s ) is a dimensionless function of the fastness parameter cu s , d is the distance to the 
system, and F is its X-ray flux. Determining the magnetic field of the neutron star from 
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a spectral analysis, we can estimate the distance to the binary system from Eq. (2) to be 
d ~ 3.1±0.7 kpc, which agrees, within the error limits, with the results of optical observations 
(Janot-Pacheco et al. 1987). 

Vela X-l 

During almost the entire period of its INTEGRAL observations, the X-ray pulsar Vela 
X-l exhibited an intense outburst activity (see, e.g., Krivonos et al. 2003) accompanied by 
significant changes in its pulse period (Fig. 1). At the same time, the shape of the source 
spectrum was almost constant, as were the positions of the two cyclotron line components 
at energies ~ 24 and ~ 50 keV (Filippova et al. 2005). Note that the absorption line at 
~ 50 keV had long believed to be the fundamental harmonic of the cyclotron line and only 
the RXTE (Kreykenbohm et al. 2002) and INTEGRAL (Filippova et al. 2005) observations 
revealed a weaker lower harmonic at energies 24-26 keV. 

The source pulse profile has a fairly complex morphology and depends strongly on the 
energy band: at low energies, the profile is an asymmetric five-peaked one with its gradual 
simplification to a double-peaked sinusoidal profile above ~ 10 keV (see, e.g., Kreykenbohm 
et al. (2002) and references therein). The evolution of the pulse profile in hard (> 20 keV) 
energy channels is shown in Fig. 8. Emission pulsations from the source are detected up to 
~ 100 keV with a double-peaked profile: the first peak is triangular in shape with a gradual 
rise and a sharp decay; the second peak is almost symmetric relative to its maximum. We 
see from the intensity map that the relative intensity of the first peak increases with energy, 
having a maximum near the first cyclotron line harmonic, while its maximum approaches 
that of the first peak in phase (Fig. 8). PF for Vela X-l is rather high compared to other 
pulsars, ~ 60% in the 20-30 keV energy band, and gradually increases with energy, having 
a local maximum near 50 keV, where the second cyclotron line harmonic is observed in the 
source spectrum (Fig. 8). 

Cen X-3 

Although the X-ray pulsar Cen X-3 is believed to be a persistent source, its intensity 
can change dramatically (by more than a factor of 10) both during the orbital cycle and 
on long time scales. These changes can be accompanied by significant changes in the pulse 
period (Fig. 9). Nagase et al. (1992) detected a possible cyclotron absorption line in the 
source spectrum at energy ~ 30 keV, which was confirmed by Santangelo et al. (1998). 
Subsequently, Burderi et al. (2000) showed that the position of the cyclotron line depends 
strongly on the pulse phase and can change with the range from ~ 28 to ~ 36 keV. Despite 
dramatic changes in the intensity of Cen X-3 during its INTEGRAL observations, we found 
no significant changes in the position of the cyclotron line with a mean energy of ~ 31 ± 1 
keV associated with them. 

The source pulse profile depends on the orbital phase and energy band: in a high state, 
the profile is double-peaked at low energies; the relative intensity of the peaks changes with 
increasing energy and the main peak disappears almost completely at high energies (above 
~ 20 keV) (Nagase et al. 1992). The X-ray pulsar 4U 0115+63 exhibits a similar behavior 
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(Tsygankov et al. 2007). PF for the source is rather high (~ 60% at 20-30 keV) and increases 
monotonically with energy without any features (Fig. 9). 

GX 301-2 

The pulse profile for the X-ray pulsar GX 301-2 has been adequately investigated over a 
wide energy range (see, e.g., Borkus et al. 1998; Tsygankov et al. 2004). It has a double- 
peaked sinusoidal shape with a constant ratio of the peak intensities at different energies. 
Analysis of the behavior of the cyclotron line revealed a strong variability of its position 
in the pulsar spectrum, depending on the pulse phase and intrinsic source luminosity (La 
Barbera et al. 2005). Thus, for example, the cyclotron line energy in a state with a low 
luminosity (L x ~ 0.8 x 10 37 erg/s) was ~ 45 keV, while it was ~ 53 keV at maximum 
luminosity [L x ~ 2 x 10 37 erg/s). Similar changes in the cyclotron energy with luminosity 
were also revealed by the INTEGRAL observatory (Filippova et al. 2005). 

Despite significant changes in the flux from the source, especially near the perias- 
tron and apoastron, its INTEGRAL pulse profile remains fairly stable (see Fig. 10 and 
http://hea.iki.rssi.ru/integral/pulsars), while PF increases with energy almost linearly. Be- 
low, we will show that, in this case, it depends significantly on the source luminosity. 

OAO 1657-415 

Despite the fairly long INTEGRAL exposure time for this pulsar and the hard spectrum 
of this object, as yet no cyclotron absorption line has been reliably detected in its spectrum 
(Filippova et al. 2005; Barnstedt et al. 2008). This may be indicative of a fairly strong 
magnetic field near the neutron star surface, just as is observed with INTEGRAL in other 
pulsars (see, e.g., Tsygankov and Lutovinov 2005a, 2005b). Note also that the pulsar exhib- 
ited a dramatic variability, both aperiodic one and one related to the orbital cycle, during 
its INTEGRAL observations. However, the pulse period changes revealed no correlation 
between them and the changes in the source flux (Fig. 11). 

The pulse profile is in the shape of a wide asymmetric peak with two maxima in the 
standard X-ray energy band whose peculiarity is a dip at phase 0.4-0.5 (see, e.g., Lutovinov 
et al. 1994). The amplitude of the second maximum decreases as we pass to the harder 
energies, the dip is essentially "smeared", and the pulse profile turns into an asymmetric 
one with a maximum near phase 0.25 (Fig. 11); PF increases with energy monotonically, 
without any features (see Fig. 11 and Barnstedt et al. (2008)). 

Her X-l 

The X-ray pulsar Her X-l is the first source with a cyclotron resonance scattering feature 
detected in its spectrum (Truemper et al. 1978). Over several decades, the properties of 
this source have been studied by many authors based on data from various X-ray observa- 
tories. However, only in recent years has serious progress been achieved in this question: in 
particular, Staubert et al. (2007) showed that, in contrast to bright pulsars, the energy of 
the fundamental cyclotron line harmonic in Her X-l is directly proportional to the source 
luminosity, while using data from the Suzaku observatory, Enoto et al. (2008) have recently 
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confirmed the presence of a second cyclotron line harmonic in the spectrum at ~ 73 keV, 
the possible existence of which was previously claimed by Di Salvo et al. (2004) based on 
BeppoSAX data. This discovery makes this object the sixth in the list of pulsars with the 
detected second cyclotron frequency harmonic. 

The pulsar Her X-l is one of the few X-ray sources with variability simultaneously on 
three different time scales: the pulse period, the orbital period, and the super-orbital period 
related to accretion disk precession. The source pulse profile depends strongly on both energy 
(Fig. 12) and phase of the super-orbital period (http://hea.iki.rssi.ru/integral/pulsars); it 
was analyzed in detail (including phase-resolved spectroscopy) on the basis of INTEGRAL 
data (Klochkov et al. 2008b). Here, we will note the behavior of PF with energy: apart from 
its increase with energy, a local feature at energies ~ 30 — 45 keV (Fig. 12) located near the 
cyclotron absorption line is present on the plot. This feature can be interpreted either as a 
local maximum shifted relative to the fundamental harmonic toward the lower energies or 
as a local PF minimum shifted toward the higher energies. 

GX 1+4 

Ferrigno et al. (2007) claimed the possible detection of a cyclotron feature in the source 
spectrum based on INTEGRAL data at ~ 34 keV. Previously, the magnetic field was es- 
timated only by indirect methods (see, e.g., Cui 1997); its strength turned out to be an 
order of magnitude larger than the above value. Note that based on first 2.5 years of the 
INTEGRAL observations, Filippova et al. (2005) detected no cyclotron feature in the source 
spectrum. 

As regards the changes in the pulse profiles, many authors (see, e.g., Greenhill et al. 
(1998) and references therein) showed that the profile clearly depends not only on the source 
intensity but also on the spin-up. According to IBIS data, the source pulse profile is single- 
peaked with a slight tendency for the separation into two sub-peaks with increasing energy 
(Figs. 13 and 14). In some cases (Fig. 13), the separation into two sub-peaks is observed 
even in the softest available energy channel. A detailed analysis of the double-peaked pulse 
profile structure based on INTEGRAL data can be found in Ferrigno et al. (2007). PF is 
practically independent of the energy and source luminosity; note only its local increase near 
the energy 45 keV (Fig. 13) in March 2003. 

Throughout the history of its studies, the pulsar GX1+4 had unique spin-up and spin- 
down characteristics. During its INTEGRAL observations, the pulsar spun down almost 
uniformly approximately until May 2004, following which the pulse period remained rela- 
tively constant. At the spin-down phase, the period increased with a mean rate of ~1.8 
s/yr. 

EXO 2030+375 

Like the sources 4U 0115+63 and V 0332+53, the X-ray pulsar EXO 2030+375 is a 
member of a binary system with a Be star and manifests itself during outbursts; in contrast 
to the two previous objects, the outburst observed from EXO 2030+375, can be both normal, 
with an intensity ~ 50 — 150 mCrab, and giant, when the flux from the source increases almost 
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to 1 Crab. During one of such giant outbursts (~ MJD53900 - 54100), the INTEGRAL 
observatory detected a significant spin-up of the neutron star (Fig. 1), which, in addition, 
correlated linearly with the rise in the source luminosity (see also Klochkov et al. 2007). 

EXO 2030+375 is one of the first X-ray pulsars to be systematically studied in the area 
of the dependence of the pulse profile on the intrinsic luminosity (Parmar et al. 1989). The 
authors found that the profile changes both with luminosity and with energy in the standard 
X-ray energy band. Significant changes in the pulse profile with source luminosity are also 
observed in the hard energy band (Figs. 15 and 16), with the energy dependences of PF 
being similar in both cases. For a bright state during an outburst, this dependence is given 
in Klochkov et al. (2008a). These authors also detected a possible cyclotron absorption 
line in the source spectrum. This feature manifests itself only in a narrow range of pulse 
phases at energy ~ 63 keV. Unfortunately, such a high energy does not allow the features in 
the dependence of PF on energy at the cyclotron frequency observed in other pulsars to be 
detected at a statistically significant level. 



CONCLUSIONS 



We analyzed in detail the timing characteristics of the hard (> 20 keV) X-ray emis- 
sion from ten bright X-ray pulsars observed by the INTEGRAL observatory. In partic- 
ular, we have investigated the dependence of PF on the luminosity and energy band for 
the first time and showed that PF increases with energy for all pulsars; in many cases, 
this increase is not monotonic, exhibiting local features near the cyclotron line harmonics, 
which is most likely attributable to the effect of resonance absorption. In addition, we com- 
piled a catalog of pulse profiles in five hard energy channels (20-30, 30-40, 40-50, 50-70, 
and 70-100 keV) for all the observed states of the pulsars under study and their luminos- 
ity levels (http://hea.iki.rssi.ru/integral/pulsars). A characteristic feature of most of the 
double-peaked pulse profiles is a decrease in the relative intensity of one of the peaks with 
increasing energy and decreasing source luminosity. Such a behavior and the increase in 
PF with energy noted above can be qualitatively explained in terms of a simple geometric 
model proposed by Tsygankov et al. (2007) and Lutovinov and Tsygankov (2008). In this 
model, accretion columns whose height depends on the pulsar luminosity are formed at the 
poles of a neutron star with a dipole magnetic field; the rotation axis of the neutron star is 
inclined with respect to the magnetic dipole axis and the observer in such a way that the 
latter can see the accretion column over its entire height at one of the poles and only its 
upper part at the other pole (see Fig. 17). Since the temperature of the accretion column 
increases toward its base (Basko and Sunyaev 1976), high-energy photons are emitted from 
regions close to the neutron star surface, while "soft" photons are formed in the upper part 
of the column. Thus, an observer will see two formation regions of soft emission and only 
one formation region of hard emission, which will correspond to a change of the pulse profile 
from double-peaked to single-peaked with increasing energy. As the luminosity decreases, 
the column height also decreases; the visible part of the second column decreases and it can 
disappear altogether from the observer's field of view at some time, which also leads to the 
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transition from a double-peaked pulse profile to a single-peaked one. The observed increase 
in PF with energy can also be naturally explained in terms of this model: at higher energies, 
the contrast between the minimum and maximum visible surfaces of the accretion columns 
is highest. Note that a similar model was proposed by David et al. (1998) to explain the 
dependence of the pulse profile for the X-ray pulsar GX1+4 on its intensity. 

It should be noted that the simple model described above gives only qualitative (in the first 
approximation) explanations for some of the observed effects. During further (quantitative) 
modeling, the temperature distribution along the column, the beam pattern of emission, its 
dependence on energy and luminosity, the bending of light rays in strong gravitational fields 
near the neutron star surface, etc. should be included in the model. 

The increase in PF with energy may also be related to the characteristic shape of the 
spectra for X-ray pulsars with an exponential cutoff at high energies. The local PF features 
near the cyclotron line harmonics may be associated with changes in the opacity of matter 
at these energies and its strong dependence on the viewing angle. In that case, the cyclotron 
line parameters (energy, optical depth, equivalent width, etc.) will change significantly over 
one rotation (pulse) of the neutron star (Tsygankov et al. 2009, prepared for publication). 

The Dependence of PF on Source Intensity 

As we showed above, the pulse profiles for X-ray pulsars often depend on their intensity, 
which, in turn, can change PF. To test and illustrate this assumption, we constructed the 
dependence of PF in the 25-45 keV energy band on the source intensity in the same band 
for the entire set of available data (in contrast to Fig. 1, the results of this analysis are 
also presented for the pulsar Her X-l, which was observed in different states, and are not 
presented for the pulsar A0535+63, for which the dynamic range of observed INTEGRAL 
fluxes was too small for any changes in PF to be observed). 

X-ray pulsars for which PF decreases with source luminosity can be separated from the 
set of results presented in Fig. 18: GX 301-2 and OA01657-415. Such a behavior of PF 
can be explained in terms of the model described above. Pulsars that exhibit a fairly large 
spread in PF at close luminosities but show no general trend in the wide dynamic range 
of observed fluxes represent another group: Vela X-l, Cen X-3, and GX 1+4. In this case, 
the source luminosity is most likely insufficient for the formation of an accretion column 
and the emission is formed near the neutron star surface (the latter is indirectly confirmed 
by the absence of any change in the position of the cyclotron line with luminosity for Cen 
X-3 and Vela X-l); the observed changes in PF may be related to local inhomogeneities of 
the stellar wind or accretion flows. For the X-ray pulsar Her X-l, the sharp decrease in PF 
observed during a low on state may stem from the fact that at this moment the emission 
is partially scattered in the corona at the edge of a precessing accretion disk (Jones and 
Forman 1976). The PF dependence obtained for the pulsar V 0332+53 is most difficult to 
explain: PF decreases with increasing luminosity at low luminosities, as is expected from 
the model described above, but it increases with luminosity at high luminosities. 
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Fig. 4: Pulse profile for the X-ray pulsar 4U 0115+63 in the (a) 20-30, 
(b) 30-40, (c) 40-50, (d) 50-70, and (e) 70-100 keV energy bands from 
IBIS/INTEGRAL data, orbit 238 (see Table 1). (f) The corresponding rela- 
tive intensity map in "energy — pulse phase" coordinates and (g) the energy 
dependence of PF. 
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Fig. 5: Same as Fig. 4 for the source V0332+53, orbit 273, a high state. 
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Fig. 6: Same as Fig. 4 for the source V0332+53, orbit 284, a low state. 
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Fig. 7: Same as Fig. 4 for the source A0535+262, orbit 352, averaged. 
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Fig. 8: Same as Fig. 4 for the source Vela X-l, orbits 373-383, averaged. 
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Fig. 9: Same as Fig. 4 for the source Cen X-3, orbit 192, averaged. 
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Fig. 10: Same as Fig. 4 for the source GX 301-2, orbit 323. 
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Fig. 11: Same as Fig. 4 for the source OAO 1657-415, orbit 364. 
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Fig. 12: Same as Fig. 4 for the source Her X-l, orbit 339, a high state. 
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Fig. 13: Same as Fig. 4 for the source GX 1+4, orbit 46. 
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Fig. 14: Same as Fig. 4 for the source GX 1+4, orbits 235-237. 
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Fig. 15: Same as Fig. 4 for the source EXO 2030+375, orbit 190, a low state. 




Fig. 16: Same as Fig. 4 for the source EXO 2030+375, orbit 462, a high state. 
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Fig. 17: Geometry of the radiation model for accretion columns of an X-ray pul- 
sar (darker regions correspond to higher temperatures) and relative positions of the 
observer, the rotation axis, and the magnetic axis (for more detail, see the text). 
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Fig. 18: PF for X-ray pulsars in the 25-45 keV energy band versus flux recorded from 
them in this energy band. 



